function path=ouprocess(alpha,m,y0,beta,N)
%%%pg 267 of shrevebook
%m=1; alpha=0; y0=0; beta=1;
time=linspace(0,1/2,N);
deltat=1/2/(N-1);
path=zeros(size(time));
path(1)=y0;
% hold on
% for pathnum=1:100
for ii=2:N
    path(ii)=path(ii-1)+alpha*(m-path(ii-1))*deltat+beta*randn(1)*sqrt(deltat);
end
% plot(time,path);
% end
